Quantum antagonism of ferromagnetic order 
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We study the interplay of superfluidity, glassy and magnetic orders of hardcore bosons with random frus- 
trating interactions. In the classical limit, this model reduces to a ±J Edwards- Anderson Ising model with 
concentration p of the ferromagnetic bonds, which hosts a glassy-ferromagnetic transition at a critical concen- 
tration pf ~ 0.77 on a 3D cubic lattice. Our quantum Monte Carlo simulation results using the worm algorithm 
show that quantum fluctuations stabilize the coexistence of superfluidity and glassy order ( "superglass"), and 
shift the (super)glassy-ferromagnetic transition to pc > pf. In contrast, antiferromagnetic order coexists with 
superfluidity to form a supersolid, and the transition to the glassy phase occurs at a higher p. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
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Introduction- Spin glasses are frustrated magnetic sys- 
tems with quenched disorder, hosting glassy phases with ex- 
tremely slow dynamics 111]. Traditionally, the simplest model 
that exhibits spin glass (SG) behavior is the ±,/ Ising model, 
or Edwards-Anderson (EA) modelf^l, in which Ising spins 
interact via randomly distributed nearest-neighbor (NN) fer- 
romagnetic (FM) and antiferromagnetic (AFM) bonds with 
probability p and 1 — p respectively. This model has been 
studied extensively with Monte Carlo simulations and ex- 
hibits a finite temperature glass transition in 3D fslQl- It is 
shown that this model has a second order transition from a 
SG to a FM (AFM) phase at a critical concentration pf ^ 
0.77 (1 — pf ^ 0.23) as one increases (decreases) the frac- 
tion of FM bonds, and there exists reentrance of the spin glass 
phase in the temperature-disorder phase diagramlSHTI]. One 
natural question to ask is how the transition can be changed 
by quantum fluctuations ifs Uloll . In particular, understanding 
the behavior of granular superfluidity in a frozen amorphous 
structure may give us hints on the microscopic mechanism of 
the formation of supersolids (SS) ll ill . 

The observation of supersolidity in solid Helium 4 (He"*) 
I Uj] has spurred immense interests on the connection between 
superfluidity and disorder in bosonic systems. Strong experi- 
mental evidences suggest that disorder may play a role in how 
the supersolid formsll2l ll3ll . A recent experiment observed 
ultraslow dynamics in solid He^ ll4ll . suggesti ng a gla ssy type 
of supersohd, or a "superglass" (SuG) ISl [lOl [iSl ll6(] . Quan- 
tum Monte Carlo (QMC) studies on the extended hardcore 
Bose-Hubbard model with random frustrating interactions on 
a 3D cubic lattice indicate that glassiness can coexist with su- 
perfluidity, and quantum fluctuations and random frustration 
are both crucial in stabilizing the superglass state|[?|. 

While the phase transitions into superfluidity or mag- 
netic/glassy states have been studied in some detail, the na- 
ture of the various transitions into a SS/SuG are yet unknown. 
In this Letter, we study the extended hardcore Bose-Hubbard 
model with random frustrating interactions model, in order 
to characterize transitions by tuning the amount of disorder 
present in the interactions, which is potentially relevant to 



the behavior noted in He^. We demonstrate that the pres- 
ence of quantum exchange terms greatly increases the tem- 
perature dependence of the SuG-FM transition by pushing the 
low temperature phase boundary to a higher critical concen- 
tration than p^'. The SuG-SS phase boundary is also drawn to 
a higher critical concentration, and an asymmetry for SuG-SS 
and SuG-FM transitions arises. 

Model- The Hamiltonian for hardcore boson with a ran- 
dom NN interaction is 

n = -Y,V,,{n,~\/2){n,-l/2)-tY^{b\b,+hb]), (1) 

where indicates the NN lattice sites, ni is the number 
operator for hard-core bosons at lattice site z, t is the hopping 
parameter and Vij are interactions with a bimodal distribution 
given as 



p{V,,) = p6{V,, -V) + {\- p)5{V,, + V). 



(2) 



This model can be readily mapped into the standard XXZ 
model: 



H — ^ JzS^Sj — -Jxy Sj' + S^Sj 



(3) 



Vij. This model reduces to the 



where J^y = 2t and Jz ~ 
classical EA model if Jxy = 0. 

There exist both diagonal and off-diagonal long-range or- 
ders in this model. Possible diagonal long-range orders are 
FM, AFM and SG, with corresponding order parameters: 
magnetization (FM), staggered magnetization (AFM), and the 
Edwards-Anderson order parameter (SG). The EA order pa- 
rameter is defined as 



(4) 



where (• • • ) denotes a thermal average and [• • • ]av an aver- 
age over disorder realizations. As the EA order parameter 
will also capture FM and AFM ordering, one must look for a 
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TABLE I. Parameters of the simulations 
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FIG. 1 . (Color online) /? vs p phase diagram. The dashed lines indi- 
cate the position of the classical transitions SG-(A)FM, and the solid 
lines with arrows show the trend of the phase boundary shifts we find 
in the quantum model. The approximate superfluid phase boundary 
(thick solid arc) is included for reference. 



non-zero while the other order parameters remain zero in 
order to identify an SG phase. 

To determine the phase transition point, we look at the 
Binder cumulants llTIl for the order parameter, which should 
cross at the transition point for different system sizes. The 
Binder cumulants of the magnetization m = -i- Sf , and 
staggered magnetization ms = j^Tlii^^y 
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which approaches one in the FM/AFM phase and goes to zero 
otherwise. We measure the superfluid density through wind- 
ing number fluctuations 111 811 . 
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where Wk is the winding number along the k direction. In 
our simulation, we identify the SS phase by the coexistence 
of AFM and SF orders, and the SuG phase by the coexistence 
of SG and SF orders. 

Method- We use the worm algorithm QMC Il9ll to solve 
this model and simulate many different realizations of the 
disorder-sets of interactions V^^-for a given parameter set. In 
the following, we choose V/t ~ A such that the superglass 
phase exists at p = O.5|0|. The average over these realiza- 
tions gives us an estimate of the thermodynamic limits of our 
observables. Each realization is simulated independently and 
equilibrium is determined by reaching measurements that sta- 
bilize within our error bars. Typically, for the hardest cases 
this requires on the order of 10^ Monte Carlo steps, where one 
Monte Carlo step is defined as Nsites completed worm loop 
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updates. Equilibration times vary across phases: worm update 
loops tend to be long while in SF phases, and under increas- 
ing ferromagnetic order the worms have a low probability of 
hopping producing short loops. After equilibration, the final 
measurements are taken for ~ 1 /5 of the required equilibra- 
tion steps. Since most of the measurement error arises from 
the disorder average, it is inefficient to spend large amounts 
of computer time getting better statistics on each individual 
sample. 

Results- The classical 3D EA model undergoes a T = 
SG-FM transition at p'f ^ 0.770, with very little temperature 
dependence at finite T. Also, the classical phase diagram is 
symmetric about p = 0.5 if one identifies the ferromagnetic 
with the antiferromagnetic regime. Figure [T] summarizes our 
QMC results in a temperature-disorder {T-p) phase diagram. 
The SuG-FM transition appears to be strongly temperature de- 
pendent, deviating from the classical behavior Likely this is 
due to the encroaching superfluid order and subsequent com- 
petition. Meanwhile, the AFM state allows for superfluidity, 
creating a region of supersolidity within the arm of the su- 
perfluid transition line, which may be why the AFM state is 
more stable in the quantum model. Furthermore, the SuG-SS 
line pushes deeper into the SG phase, suggesting interesting 
correlation between the AFM and SF orders. 

The results for our analysis of the magnetic ordering are 
presented in Fig. |2] For an even mix of bonds p = 0.5, the 
transition to superfluidity occurs near /3c ^ 1. Given that we 
expect the superfluidity to be killed off by ferromagnetism, it 
is reasonable to guess that j3c should increase as p increases. 
Thus, for our data in the vicinity of Pc\ at /3 = 1, we don't 
expect much interference from quantum mechanical effects. 
Data at /? = 1.5, 2 (not shown) are similar, though with finite- 
size effects (FSE) starting to increase in magnitude. 

At /? = 3 in Fig. |2] (a), it is also clear there would not 
be a transition for p < 0.80. This is markedly different 
than the classical results. The classical transition line pf as 
a function of /3 is only weakly temperature dependent, shift- 
ing less than a percent between the multicritical point and zero 
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FIG. 2. (Color online) (a) Finite size scaling of Qm at /3 = 1, 3 to determine the SG-FM transition. Data for /? = 3 are shifted vertically for 
clarity. 13 — 1 agrees with the classical model while /3 = 3 suggests the transition is pushed to higher p. (b) Data for ps showing the transition 
out of superfluidity appearing in the same vicinity. 
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FIG. 3. (Color online) Finite size scaling of the staggered magnetiza- 
tion cumulant Qsm- Data for /? = 2 are shifted vertically for clarity. 
Results at /? = 1 suggest the SuG-SS transition has been drawn into 
the classical SG region. 



temperature Here, we show at least a four percent change 
at /? = 3 compared to pf. We have studied the /? = 3 su- 
perfluid transition as well, and it appears in the same vicinity 
[Fig. |2] (fe)]. We mention that the scaling form ^ 
deviates from the expected scaling form, and this is addressed 
separately ll23l . 

The specific mechanism behind the increased temperature 
dependence of the SuG-FM transition should be rooted in how 
superfluidity favors the spin glass phase over the ferromag- 
netic phase. We hypothesize that, first, as the exchange terms 
cause the system to destroy ferromagnetic order to allow hop- 
ping of particles, clusters of FM-ordered spins are suppressed. 
Second, glassy clusters more readily allow superfluid fluctua- 
tions amongst them, perhaps following lines of sites that are 




FIG. 4. A sample arrangement one might see in a state near the 
SuG-SS line. The hop move indicated is energetically neutral with 
respect to the bonds and destroys the local AFM order. However, the 
new state allows fewer fluctuations so is not favored with respect to 
entropy. 

weakly constrained due to frustration. Finally, this results in 
the transition line moving to allow for a broader spin glass 
(thus superglass) phase. 

We also looked at the behavior on the other side, p < 0.5, 
which is different due to the asymmetry of the ground state 
quantum fluctuations-spin exchange is forbidden out of a pure 
FM ground state, while possible from a pure AFM ground 
state. Notably, we find superfluid transitions at higher temper- 
atures: /?c(l— p) < Pdp)- Fig. [3] shows data for the Binder 
cumulant using the staggered magnetization, gsm at /3 = 1 
and 2. The staggered magnetization data suffer from larger 
FSE and subsequently require larger lattices to achieve com- 
parable precision, restricting the current study from exploring 
lower temperatures. 

The position of the crossing at /3 = 1 lies around pc ^ 0.25 
indicating that the SuG-SS line is noticeably shifted from clas- 
sical behavior The data at lower T are less conclusive due to 
larger sizes being out of reach at present, but suggest an even 
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larger shift. Above the SF transition at /3 = 0.7(T = 1.43), 
the AFM-SG transition agrees with the classical transition 
point p = 1 — pf = 0.23. Thus, the stronger the quantum 
mechanical nature becomes-the deeper into the superfluid be- 
havior we go-the stronger the AFM order appears. This is at 
first counterintuitive, but we can provide a simple, consistent 
picture for both this and the SuG-FM behavior When an FM 
bond is satisfied, i.e. two neighboring sites are in the same 
state, this forbids hopping along that bond. Conversely, when 
an AFM bond is satisfied, the two sites are in an opposite state 
and will allow the system to gain kinetic energy through hops. 
As this is energetically favorable, we can consider the AFM 
bonds to have a larger effective bond strength. Consequently, 
the fraction of frustrating FM bonds it would require to de- 
stroy the AFM order should rise, and we see the present shift 
in phase boundary. 

Speaking more generally, we can consider the change in 
the entropy of the system, with respect to quantum fluctua- 
tions, when moving towards higher AFM or FM order Under 
pure FM order no hopping is allowed, while under pure AFM 
order each site has the ability to engage in virtual hopping. 
Using a 2D lattice as an example, Fig.|4]shows a typical local 
snapshot of a system with AFM order being frustrated by FM 
bonds. The indicated move would go against three NN bonds 
while satisfying three others, making the diagonal component 
energetically neutral, and this would turn the local AFM or- 
der into FM order However, after this move the system is 
now considerably more constrained regarding the number of 
hoppings available. Originally, the two interior sites could be 
involved with seven different hop moves, but these sites now 
have a single move: to return back to the original state. In this 
way, the AFM state is favored as it allows for more gain in 
kinetic energy. 

Summary- We have presented our QMC results of a 
bosonic model exhibiting a superglass phase in order to study 
the phase transitions achieved by directly tuning the level of 
disorder present. These results indicate that the addition of the 
exchange terms act to stabilize the classical spin glass phase 
against the formation of ferromagnetic clusters which impede 
hopping. In addition, the favoring of AFM bonds to FM bonds 
leads to a shift in the SuG-SS phase boundary. An interesting 
issue not addressed here is whether there is alwa ys a SG phase 
between the FM and SuG phases. Recent work i2lll on disor- 
dered Bose systems suggests that this may be true. Unfortu- 
nately, our current precision is not high enough to be conclu- 



sive. Future work is required to clarify this issue. 

We thank Anders Sandvik and Markus Miiller for helpful 
discussions. We are grateful to National Center for High- 
Performance Computing, and Computer and Information Net- 
working Center, NTU for the support of high-performance 
computing facilities. This work was partly supported by the 
NSC in Taiwan through Grants No. 100-21 12-M-002 -013 - 
MY3, 100-2120-M-002-00 (Y.J.K.), and by NTU Grant num- 
bers 10R80909-4 (Y.J.K.). Travel support from National Cen- 
ter for Theoretical Sciences is also acknowledged. 



* |yjkao@phys.ntu.edu.tw| 
[1] K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 (1986). 
[2] S. F. Edwards and P W. Anderson, J. Phys. F: Met. Phys. 5, 965 
(1975). 

[3] N. Kawashima and A. P Young, Phys. Rev. B 53, R484 (1996). 
[4] H. G. Ballesteros etai, Phys. Rev. B 62, 14237 (2000). 
[5] A. K. Hartmann, Phys. Rev. B 59, 3617 (1999). 
[6] M. Hasenbusch, A. Pelissetto, and E. Vicari, Phys. Rev. B 78, 
214205 (2008). 

[7] G. Ceccarelli, A. Pelissetto, and E. Vicari, 

[Phys. Rev. B 84, 134202 (201 1)| 
[8] T. E. Markland, J. A. Morr one, B. J. Berne, K. Miya zaki, E. Ra- 

bani, and D. R. Reichman, [Nat Phys 7, 134 (201 1)| 
[9] K.-M. e. a. Tarn, Phys. Rev. Lett. 104, 215301 (2010). 
[10] X. Yu and M . Mueller,_ ArXiv e-prints (2011), 

larXiv: 1 1 1 1.5956 [ cond-mat.d is-nn]| 
[11] E. Kim and M.H. W.Chan, Nature 427, 225 (2004); 

E. Kim and M. H. W. Chan, Science 305, 1941 (2004)] 

|http://www.sciencemag.org/content/305/5692/1941.full.p df 
[12] A. S. C. Rittner and J. D. Reppy, 

[Phys. Re v. Lett. 98, 175302 (2007) 
[13] A. S. C. Rittner and J. D. Reppy, 

[Phys. Rev. Lett 97, 165301 (2 006) 
[14] B. Hunt, E. Pratee, V. Gadagkar, M. Yamashita, A. Balatsky, 

and J. C. Davis, Science 324, 632 (2009). 
[15] G. Biroli, C. Chamon, and F. Zamponi, 

IPhys. Rev. B 78, 224306 (2008) 

[16] J. Wu and P Phillips, Phys. Rev. B 78, 14515 (20087[ 
[17] K. Binder, Z. Phys. B 43, 119 (1981). 

[18] E. L. Pollock an d D. M. Ceperley, 

[Phys. Rev. B 36, 8343 (1987)| 
[19] N. Prokof ev, B. Svistunov, and I. Tupit- 

syn, Physics Letters A 238, 253 (1998) , 

[Journal of Experi mental and Theoretical Physics 87, 310 (1998)^ 
10.1134/1.558661^ ' 
[20] D. A. Larson and Y.-J. Kao, Unpublished (2012). 
[21] L. Pollet, N. V. Prokof ev, B. V. Svistunov, and M. Troyer, 
|Phys. Rev. Lett. 103, 140402 (2009)] 



